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Abstract: 

We consider the inhomogeneous decay of an unstable D-brane. The dynamical 
equation that describes this process (in light-cone time) is a variant of the non- 
linear reaction-diffusion equation that first made its appearance in the pioneering 
work of (Luther and) Fisher and appears in a variety of natural phenomena. 
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1 Introduction 



D-branes play an important role in the dynamics of and duality relations between different 
string theories. These are defined by boundary conditions on open strings, but at the same 
time are typically stringy non-perturbative excitations of closed strings. While some D- 
branes carry conserved charges and are stable, others are unstable. To be precise, all D- 
branes of different dimensions of the bosonic string theory and more than half of those of 
the type I and II superstring theories are unstable. Moreover, configurations of more than 
one kind of D-branes could be unstable even when the constituents are stable individually; 
for example, a brane-anti-brane pair. A detailed understanding of this instability provides 
a valuable window into the behaviour of strings. 

Even though a string has an infinite number of vibrational modes, the leading effects 
due to the instability and its characteristic features are governed by the lowest mode, the 
so called tachyonic scalar field. The study of unstable branes was pioneered by Sen [T] , who 
proposed a set of precise conjectures concerning the tachyon. Since then, these conjectures 
have been checked in the open string field theory in the level truncation scheme, in various 
toy models, in the boundary string field theory, and finally, an exact solution describing 
the state at the (local) minimum of the tachyon potential has been constructed in open 
string field theory. Much less studied, however, is the time- dependent dynamical process 
of tachyon condensation [2HS]. 

Ignoring all the details and and technicalities, this is a system which has two extrema: 
an unstable maximum (the perturbative vacuum of open strings) and a stable (or locally 
stable, as in the case of bosonic string) minimum. It is expected that the system will 
make a transition from the unstable to the stable phase dynamically. This situation is 
ubiquitous not only in physics but in various other fields, among them biological and 
chemical systems. A typical equation that governs the dynamics in such cases is the 
Fisher equation: 

d t u(t, x) = Dd?u(t, x) + ru(t, x) (1 — u(t, x)) , (1.1) 

where, D and r are constants; (a more general function f(u) may be considered in the 
RHS). It admits a time-dependent solution that corresponds to a front which separates 
the two phases (an unstable one at u = and a stable phase at u = 1) and moves with 
a characteristic speed while retaining its profile. This reaction-diffusion equation and its 
travelling front solution has a long history: Although it was first written by Luther (1906) 
for a chemical system, unaware of this work, Fisher (1937) proposed this equation and 
studied its front solution to describe the spread of an advantageous mutation. A detailed 
mathematical analysis by Kolmogorov, Petrovsky and Piskunov was the first in the vast 
literature [T0"l[TT] that followed. The equation of the tachyon field on an unstable D-brane 
turns out to be a variant of this with new elements in the form of time delay and spatial 
non-locality. 
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2 Open string field theory 



Henceforth, for definiteness, we shall restrict to the unstable branes of the bosonic string 
theory. The tachyon equation can be obtained from the cubic open string field theory. 
The string field, expanded in terms of the infinite number of oscillatory modes, is 



\*) = (<f>(X) Cl + ---)\0) 



d n k 
(2ir) n 



(j>{k) e lk - x Cl + 



|0> 



where <ft denotes the scalar field corresponding to the lowest mode of the string and the 
dots denote the higher excitatations that have been omitted. The cubic action of the 
string field is of the Chern-Simons type: 
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In the above, the products are defined in the conformal field theory on the upper half- 
plane in the usual fashion. If we retain only the tachyon field (level truncation to zeroth 
order) we obtain the action [2]: 



S 



d n x 



where, K = 3v^3/4. The equation of motion for the tachyon 

a' x) = -</>(t,x) + K 3 e aD [e aD 0(t,x 



(2.1) 



(2.2) 



(where, a = a' In K) contains an infinite number of higher derivatives in the interaction 
term and is, therefore, non-local. 

When (ft depends only on time (spatially homogeneous decay): 



d 2 <p(t) 
dt 2 



4>(t) 



(2.3) 



This equation has solutiona3 that start at the maximum of the potential (at <pu = 0) 
towards the (local) minimum (at <fi = K~ 3 ~ 0.456), but it overshoots and exhibits (non- 
linear) oscillations around the minimum. At late times, these behave wildly. However, 
there is not a solution that interpolates between the unstable and the stable extrema [2]. 
On the other hand, according to the conjectures of Sen, one expects the unstable D-brane 
to decay into some configuration of closed strings, which will carry the energy (density) of 
the brane. At the tree level OSFT, however, open strings do not interact with the closed 
string modes for this to happen. 



1 Empirically, the initial value problem turns out to be well- defined with just two initial conditions, 
say, the position and the velocity j2]- For further analysis of initial conditions of equations of this type, 
see Refs. [5j[9]. 
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3 Branes in Linear Dilaton Background 



In order to circumvent this, yet not deal with the complexities of an open-closed string 
field theory, one may consider open string field theory in the presence of a closed string 
background. Perhaps the simplest of these is a linear dilaton background^] considered 
in Ref . [3] . These authors use light-cone coordinates x ± = (t ± x) / y/2, and consider the 
dilation profile T>(x) = —D + x~ = —bx~ to study the homogeneous decay of the tachyon as 
a function of light-cone time x + , which we shall call r to simplify notation. The dilaton, 
being linear along a null direction, changes the (world-sheet) conformal dimension of 
the tachyon vertex operator e %k,x from k 2 to k 2 + ibk~ (but does not alter the matter 
contribution to the central charge). Consequently, the equation of motion for the tachyon 
gets modified from Eq. (12.21) to 



a '{ b ^~ ^ r ' x ^) = ^ r ' x ^) - ^ 3 e" 2 ^ +QV i [e aV i0(r,x ± )] 2 , (3.1) 

where xi denotes the coordinates transverse to the light-cone coordinates. This is the 
'Fisher equation for the tachyon on a decaying brane'. While there are many variants of 
the Fisher equation ^Ujlll]. this particular incarnation with a nonlocal nonlinear term 
and a delayed dependence on time, is, to our knowledge, noveljf] 

The case where <fi = </>(r) depends only on time (homogeneous decay), and studied in 
Refs. PEHS], 

a'b d + 4>(r) = 0(r) - K 3 [0(r - 2ab)} 2 , (3.2) 

is a canonical example of a delayed growth model used, e.g., in population dynamics [TO] . 
While the standard logistic growth model has a simple interpolating solution, the delay 
leads to oscillations around the stable fixed point at <ps- This follows from linearizing 
Eq. (13.21) around the stable fixecQ point in terms of <p = <ps + if>- 

a'bd T ip(T) = if)(r) - ip(r - 2ab). (3.3) 

An ansatz on the form if) ~ e~ Xr leads to a transcendental equation 

a'b\ + 1 = 2e- 2ab \ (3.4) 

that does not have a real solution, but an infinite number of complex solutions (occurring 
in complex conjugate pairs). This property is characteristic of delayed differential equa- 
tions. The (oscillatory) convergence to <frs is determined by the solution with the smallest 
value of Re A. 



2 Another possibility that has cosmological implications is to couple the tachyon to gravity |12) . 
3 Similar forms of nonlocality in interactions in biological systems |13H15j was pointed out to me by 
V.M. Kenkre. 

4 The linearized equation around the unstable fixed point <fiu — 0, solved by exp(r/6), shows the 
system moving exponentially away from it, as one would expect of a tachyonic scalar field. 
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4 The Travelling Front 



Getting back to the Fisher equation ( 13.11) . let us consider only one transverse coordinate 
(denoted by y) for simplicity. We seek a travelling front solution that moves from the right 
to the left (so that at any instant of time the region to the right of the front is converging 
to the stable fixed point). First consider the equation linearized around <fru = 0. A trial 
solution of the form ~ exp (k(y — v(k)r), leads to the dispersion relation 

<Q = Uk + ±), (4.1) 



that has a minimum at f m ; n = 2/6. The wavenumber k is real for v (k) > t> m in, therefore, 
any of these would solve the linearized equation. For the standard Fisher equation, with 
a large class of nonlinear interactions, the travelling front is proven to select t> m i n among 
this PHttl] • This feature is likely to be true of Eq. (13. ip (with the additional elements of 
delay and nonlocality) because the 'initial condition' — the leading edge of the wave — is 
determined by the 'mass' of the tachyon. However, we shall not dwell on a more rigorous 
proof here. 

Rather, we look for a travelling front solution in the form of 0(r, y) = = y + v r), 
and use singular perturbation analysis [TQl[Tll[T6] to determine the solution $(rj). Note, 
from Eq.f l4.ll) . that v 2 b 2 > 4. In the absence of a naturally small parameter, the idea is 
to scale 7] = yfe!; by e = l/v 2 b 2 < 0.25, so that the equation takes the form: 

- -==& 2 $ = Vtf$ - ~^K 3 exp (-2ad^ + eadfj e £a9 l<& , (4.2) 

oo 

Now expand $(£,e) = $ n (£)e n as a power series in e and compare terms. The lowest 

n 71=0 

order equations that determines $o(0 

d^ - $ + K 3 e- 2a9 t [$ (0] 2 = 0, (4.3) 

is identical to the homegeneous equation (13. 2p . Therefore, the solution of Ref. [3] is a seed 
for the travelling front. The correction at 0(e), can be solved from 

<%$i - $i + 2K 3 e- 2a9 t [$ $i] = (4.4) 
<9f $ - K 3 e~ 2a9 z [4a$ 9|$ + 2a (<%$ ) 2 ] , 
after substituting $o(0 from Eq. (I4.3|) . The equation at 0(e 2 ) is 

d^ 2 - $ 2 + 2K 3 e~ 2a9 z [$ $ 2 ] = (4.5) 
<9f$i - K 3 e- 2a9 z [4a 2 $ <9|$o + 8a 2 <%$ <9f $ + 6a 2 (<9f$ ) 2 

4a(9|$o$i + 4a<9 ? $ 9 c $i + 4a$ <9 2 $i + $?] , 



5 We shall set a 1 = 1 from now on. 
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by substituting for $o(£) and 3>i(£) from Eqs. fl4.3p and (14. 4p . respectively, and so on 
recursively. 

The functions $ n for the usual Fisher equation can be found analytically. For the 
tachyon, the propagating front can be obtained by NDSolve for delayed differential equa- 
tion in mathematica. We choose an 'initial' configuration Ae^ and adjust A such that at 
£ = 0, $(0) is half-way to the stable vacuum at (fis- Due to the delay, the tachyon settles 
to the stable vacuum after damped oscillations. For the correction at 0(e), to tame any 
unnatural behaviour of the numerical solution, we choose an 'initial configuration' for $i 
that is identical to the standard Fisher equation for £ < 0. Indeed, in this region, the 
profile of the front in the two cases are rather similar. The results are displayed in Fig{TJ 

O $o + 0.25 €>! 



0.6 r 




Figure 1: On the left: The leading order solutions <3?o(£) of the ordinary Fisher equation (gray 
dotted) and the tachyon equation (|4,3p (blue). On the right: Solutions upto 0(e) — the 0(1) 
solution of the tachyon equation is shown as blue dashed line. The undulation of the tachyon 
around the stable vacuum is characteristic of the delay. 

In adapting singular perturbation theory to our problem, we have naively truncated 
the infinite number of derivatives in e a « to a small finite number. One could proceed in 
another way to avoid this problem. We notice that the interaction term involves: 

e ad !f(0 ee <5[f] = -±= r e-^f(C)dC (4.6) 

and folding by the Gaussian kernel softens the oscillations. Since, in the limit, a — > 0, the 
kernel becomes the Dirac 5-function, u(J5[$o] = G5[3?o] — is 0{e) in our case. The lowest 
order equation (I4.3P remains unchanged, but the first order correction is now determined 
by: 

- $! + 2K 3 e~ 2ad t [$0$!] = gi ($ ), (4.7) 

where, #i($ ) = df$o - K 3 e~ 2a % (2$ 3lS[$o] + f Let us find $ i in a different 

way: We begin with the standard Fisher equation, which corresponds to a = 0. In that 
case, the analogue of the LHS of Eq. (l4.7p (and indeed all the equations at higher orders 
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in e) are simplified by the observation that 1 — 2K 3 Qq = ^o/^o- This helps to reduce the 
problem of finding $ n >o to one of quadrature. Formally, this is still true as an operator 
equation, for a differentiation of Eq.gSJ yields [l - 2K 3 ( e - 2Q ^$ ) e~ 2ad t] <9 C $ = <9f $ - 
One can now integrate in mathematica to find $1. The constant of integration is chosen 
so that the correction vanishes at £ = [ID] . The solution upto 0(e) obtained this way 
from Eqs. fl4.4p and (14. T[) are shown in Fig[2] Even though the error seems greater, either 
due to the operator identity not being too accuarte, or due to numerical integration, the 
singular perturbation method using the integral transform (j4.7p (instead of truncation) 
clearly shows improved behaviour. 




Figure 2: On the left: Solutions at 0(e) obtained by integration of Eqs. (|4.4j) (grey dot-dashed) 
and (|4.7j) (blue). On the right: Solutions upto 0(e). The 0(1) solution is shown as gray dotted 
line. 



5 Conclusions 

We close with a couple of comments. The decay of an unstable brane will be triggered by 
the tachyon moving away from the maximum of the potential in a finite region of space. In 
the one dimensional case we have studied, this will lead to two fronts, one moving to the 
left and the other to the right. In higher dimensions, in the spherically symmetric case, 
the Laplacian in Eq. fl3.ip in the radial variable r does not give a Fisher type equation, 
but will asymptote to one for large r [10]. 

In summary, the dynamical equation of a tachyon on an unstable D-brane is a Fisher 
type react ion- diffusion equation, in which the interaction is smeared by a Gaussian kernel 
and is also delayed. It will be interesting to see if these additional features are useful 
elsewhere; for example, its effect on pattern formation in biological systems [15] may be 
worth studying. As for the decaying brane, extension of the travelling front to a solution 
in string field theory as well as its stability are among the open problems. 
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